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Using the most recent data from the WMAP, ACT and SPT experiments, we update the con- 
straints on models with oscillatory features in the primordial power spectrum of scalar perturbations. 
This kind of features can appear in models of inflation where slow-roll is interrupted, like multifield 
models. We also derive constraints for the case in which, in addition to cosmic microwave obser- 
vations, we also consider the data on the spectrum of luminous red galaxies from the 7th SDSS 
catalog, and the SNIa Union Compilation 2 data. We have found that: (i) considering a model with 
features in the primordial power spectrum increases the agreement with data with the respect of 
the featureless "vanilla" ACDM model by A.y^ ~ 7; (ii) the uncertainty on the determination of the 
standard parameters is not degraded when features are included; (iii) the best fit for the features 
model locates the step in the primordial spectrum at a scale k ~ 0.005 Mpc~^, corresponding to the 
scale where the outliers in the WMAP7 data at f = 22 and £ — 40 are located.; (iv) a distinct, albeit 
less statistically significant peak is present in the likelihood at smaller scales, with a Ax^ — 3.5, 
whose presence might be related to the WMAP7 preference for a negative value of the running of 
the scalar spectral index parameter; (v) the inclusion of the LRG-7 data do not change significantly 
the best fit model, but allows to better constrain the amplitude of the oscillations. 

PACS numbers: 98.80.Cq, 98.70.Vc, 98.80.Es 



I. INTRODUCTION 

The inflationary paradigm is an integral part of the 
currently accepted concordance cosmological model, ex- 
plaining the flatness and homogeneity of the observed 
Universe, as well as providing a mechanism to produce 
the primordial curvature perturbations that eventually 
led to the formation of structures. The shape of the 
power spectrum of primordial perturbations can be con- 
strained, at least at the largest scales, using cosmic mi- 
crowave background (CMB) data. The 7-year WMAP 
data are in excellent agreement with the assumption of 
a nearly scale-invariant power spectrum of scalar per- 
turbations [11 E] ■ Such a spectrum, described by a sim- 
ple power law with spectral index ris very close to (al- 
beit different from) unity, is the one that would be pro- 
duced in the simplest inflationary scenario, that of a sin- 
gle, minimally-coupled scalar field slowly rolling down a 
smooth potential. The expectation of a power-law spec- 
trum continues to hold up against scrutiny also when 
tested against observations at scales smaller than those 
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probed by WMAP, like the small-scale CMB measure- 
ments of the Atacama Cosmology Telescope (ACT) [3HS] 
and South Pole Telescope (SPT) fS|, and the spectrum 
of luminous red galaxies [7 . Nevertheless, a scale invari- 
ant power spectrum with = 1 could be easily put in 
agreement with data in some non-minimal models, e.g. 
considering an extended reionization process [8 10], non 
standard processes during recombination like dark mat- 
ter annihilation pTMl4j . extra relativistic particles (see 
e.g. [inilH]) and so on. 

In spite of this, however, models with localized "fea- 
tures" in the primordial power spectrum provide a better 
fit to the data with respect to a smooth power- 

law spectrum. This is mainly due to the presence, in the 
WMAP temperature anisotropy spectrum, of two out- 
liers in correspondence oi £ — 22 and £ = 40. In par- 
ticular, these "glitches" are well fitted by a primordial 
power spectrum featuring oscillations localized in a suit- 
able range of wave numbers. On the other hand, it is 
worth noticing that the "glitches" could have a more con- 
ventional explanation, steaming from some still unknown 
systematics in the WMAP data. 

Features in the primordial power spectrum can be gen- 
erated following departures from slow roll, that can hap- 
pen in more general inflationary models. In particular, 
in multifield supergravity- or M-theory-inspired models 
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[231 [2S], a field coupled to the inflaton can undergo a 
symmetry-breaking phase transition and acquire a vac- 
uum expectation value. Such a phase transition corre- 
sponds to a sudden change in the inflaton effective mass 
and can be modeled as a step in the inflationary po- 
tential. The presence of the step produces, in turn, a 
burst of oscillations in the power spectrum of curvature 
perturbations [5^1 [57], localized around the scale that 
is crossing the horizon at the time the phase transition 
occurred. Departures from the standard power-law be- 
haviour can also be present in trans-planckian models 
[5SH3T] , in models with a phase of fast roll [32] , or with 
a sudden change in the speed of sound [55H55] . Simi- 
larly, in the so-called Starobinsky model [36] , a change in 
the slope of the potential causes a step in the perturba- 
tion spectrum. In addition to their effect on the power 
spectrum, these non-standard inflationary scenarios can 
also be constrained through their predicted bispectrum 

[351 [371 [3i. 

The purpose of the present work is to use current data 
to update previous constraints that have been put on the 
presence of such a step-like feature in the inflaton poten- 
tial. We improve over previous works by using a more 
complete dataset that includes the WMAP temperature 
and polarization data, the small-scale CMB data from 
ACT and SPT, and the matter power spectrum obtained 
from the Luminous Red Galaxies (LRG) sample of the 
Sloan Digital Sky Survey (SDSS) 7th data release [35]). 
The inclusion of different datasets allows us to explore 
a wider range of scales with respect to previous analy- 
ses, going from the Hubble radius down to the smallest 
linear scales, k ~ 0.1 Mpc~^. In particular, this leads 
to the clear identification of a "forbidden" range where 
oscillations are not allowed. 

The paper is organized as follows: in Section [ll] we 
briefly recall the theory concerning the evolution of infla- 
tionary perturbations in interrupted slow roll models; in 
Section III we describe the phenomenological model used 
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to describe a step in the inflationary potential, and the 
analysis method adopted in the present work; in Section 
rv| we present the results of the analysis , and in Section 
VTwe derive our conclusion. 
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to determine the background dynamics of the Hubble pa- 
rameter H and of the (unperturbed) inflaton field 0. 

In order to study the evolution of the curvature per- 
turbation 7?,, one introduces the gauge-invariant quan- 
tity [iOHi^ u = —zTZ, where z = a(j)/H and a is the 
scale factor. The Fourier modes Uk of u evolve according 
to (primes denote derivatives with respect to conformal 
time r/): 



(3) 



In the limit fc^ ^ z"/z, the solution to the above equa- 
tion should match the free-field solution — e~^'''^ / V2k. 
The evolution of z is determined directly by the solution 
of Eqs. ([T]) and (|2|, although during slow roll one can 
approximate z"/z ~ 2a^H^. At this point, it is possi- 
ble to integrate Eq. ([s]) to get Uk{r]) for free-field initial 
conditions. 

Finally, the power spectrum of the curvature pertur- 
bation P-n is related to u and z through 
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evaluated when the mode crosses the horizon. 



B. Models with interrupted slow roll 

In the following we shall consider models where slow 
roll is briefly violated. Phenomenologically, these can be 
described by adding a step feature to a V((j)) = rnP'cj)^ /2 
chaotic potential, i.e., by considering a potential of the 
form 



1/(0) = ]^m'f 
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II. INFLATIONARY PERTURBATIONS IN 
MODELS WITH INTERRUPTED SLOW ROLL 



A. Inflationary pertubations 

Let us start by briefly recalling how to compute the 
spectrum of primordial perturbations for a given infla- 
tionary potential V{4)) [26". In the following we shall 
work in reduced Planck units (c = Ti = SvrG = 1). The 
flrst step is to solve the Friedmann and Klein-Gordon 
equations (dots denote derivatives with respect to the 



where h is the value of the field where the step is located, 
c is the height of the step and d its slope. Although the 
underlying potential is taken to be the one of chaotic 
inflation, we shall see below that this form can also be 
used to describe different kinds of potential. 

A sharp step in the inflaton potential, like that de- 
scribed by Eq. ([s]), can appear for example in multi-fleld 
inflation models, following a symmetry-breaking phase 
transition undergone by another fleld coupled to the in- 
flaton. This induces a rapid variation in the inflaton 
effective mass that is reflected in the potential (in- 
deed, the potential ([s]) is of the form 1/(0) = ^ml^cj)^, 
).T^n this regard, one can think of b 
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FIG. 1. Model likelihood (upper panel) and A^^ (lower panel) 
as functions of b for the CMB dataset, obtained by maximiza- 
tion. 

as being related to the tinie when the phase transition 
occurs, c to the change in the inflaton mass, and d to the 
width of the transition. 

The spectrum of primordial perturbations resulting 
from the potential ([s]) can be calculated as outlined in 
the previous section, and is found to be essentially a 
power-law with superimposed oscillations. The oscilla- 
tions are localized only in a limited range of wavenumbers 
(centered on a value that depends on b) so that asymp- 
totically the spectrum recovers the familiar k'^'~^ form 
typical of slow-roll inflationary models. In particular, for 
a chaotic potential, the underlying power law has a spec- 
tral index ~ 0.96. 

One issue that we have left aside so far is how to relate 
the horizon size at the time the step occurs to a physical 
scale. This depends on the number of e-folds taking 
place between the time a given mode has left the horizon 
and the end of inflation. We choose Ni, = 50 for the 
pivot wavenumbcr k^, — kg — 0.05 Mpc^^. This choice 
is somewhat arbitrary; however, a different choice would 
correspond to a translation in the position of the step in 
(j) and would thus be highly degenerate with b. For this 
reason we do not treat iV^ as a free parameter, consistent 
with what has been done in previous studies [IHl US] . 

III. ANALYSIS METHOD 

In order to compare the theoretical predictions for the 
constraints on the parameters characterizing inflationary 
models with a step in the inflaton potential with observa- 
tional data, we performed a Monte Carlo Markov Chain 
analysis via the publicly available package CosmoMC [43j. 
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FIG. 2. Primordial power spectrum for an inflationary poten- 
tial of the form |5| with m — 7.5 x 10"*^. The values of the 
step parameters are : b — 14.66, logc — —2.75, logd = —1.42 
(red), b = 14.00, logc = -2.66, logd = -0.54 (magenta dot- 
ted), corresponding to the two peaks in the likelihood. For 
comparison, we also show the best-fit ACDM power spectrum 
(black dot-dashed). We note that the model with 6 = 14 re- 
sembles, in the fc— range considered, a model with a negative 
running index. 



We used a modified version of the CAMB ([S]) code in 
which we numerically solve Eqs. ([l])~(|3]) using a Bulirsch- 
Stoer algorithm in order to theoretically calculate the ini- 
tial perturbation spectrum Q, needed to compute the 
CMB anisotropics spectrum for any given values of the 
parameters describing this type of inflationary model. 
Then we compare these theoretical models with two dif- 
ferent combination of data sets. We will briefly come 
back on describing the principal characteristics of each 
of the dataset considered in this work. 

We consider chaotic inflation potentials of the type of 
Eq. [5j Following the prescription described in Sec. |II A[ 
this potential leads to a well-defined primordial pertur- 
bation spectrum V-n- The free parameters in Eq. [5] are 
then the inflaton mass m and the step parameters 6, c 
and d. In our analysis we map the mass m onto A^, i.e. 
the amplitude of the primordial spectrum at the pivot 
wavenumbcr ko = 0.05 Mpc~^, as indeed the inflaton 
mass sets the overall scale for the potential and conse- 
quently for the amplitude of the perturbations. We note 
also that the choice of the pivot wavenumbcr changes the 
relationship between the value of b and the position of os- 
cillations in A:— space; this should be taken into account 
when comparing the results of different studies. In par- 
ticular, changing feg from 0.05 to 0.002 Mpc^^ shifts b by 
~ 0.5 towards lower values. 

As previously noted in Sec. |IIB[ for the chaotic poten- 
tial of Eq. [5] the smooth power law has a fixed spectral 
index, that is Ug — 0.96. However, as noted in Ref. jl9j . 
more general forms of the potential can be phenomeno- 
logically taken into account by promoting rig back to a 
free parameter and defining a "generalized" primordial 
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FIG. 3. Upper panel: Power spectrum of CMB Anisotropies 
Power Spectrum for the best-fit ACDM model (black line) 
and two step models with 6 = 14 (red dashed) and h — 14.66 
(blue dot-dashed), corresponding to the two minima in x^, 
compared with WMAP7, ACT and SPT data. Lower panel: 
The same as above, but plotted in terms of residuals with 
respect to the ACDM best fit. 



spectrum as 



n3-0.96 



(6) 



where '/^(fc) is the spectrum induced by the chaotic po- 
tential ([fj). 

Therefore the theoretical model we are considering is 
described by the following set of parameters: 



(7) 



where = J7f,/i^ and ll>c — ^c^^ are the physical 
baryon and cold dark matter densities, 9 is the ratio be- 
tween the sound horizon and the angular diameter dis- 
tance at decoupling, r is the optical depth to reioniza- 
tion, 6, c and d are the parameters of the step-inflation 
model, As is the overall normalization of the primordial 
power spectrum (equivalent to specifying as discussed 
above), and Ug is the effective tilt. We consider purely 
adiabatic initial conditions, impose flatness and neglect 
neutrino masses, and limit our analysis to scalar pertur- 
bations. 



We consider as our basic dataset a combination of three 
different CMB datasets: WMAP7, ACT and SPT (in 
what follows we will refer to this case simply as "CMB"), 
and then we also consider an enlarged dataset, dubbed 
"CMB-hLRG", where we also add the LRG SDSS cat- 
alog 135], the Supernovae la Union Compilation 2 data 



FIG. 4. The same as Fig. [3] but zoomed in the region 700 < 
£ < 2000. Note that in order to improve clarity, the vertical 
axis shows ^*Cf instead than £{£ + l)Ci. 



06* , and impose a prior on the Hubble constant from the 
Hubble Space Telescope (HST) [47]. 

To compute the likelihood of the data we have properly 
modified the CosmoMC package in oder to make use of the 
routines supplied by the WMAP and ACT teams for their 
datasets, both publicly available from the LAMBDA 
websit^ and of the likelihood code provided by the SPT 
team |3] for the SPT dataset. 

The ACT and SPT experiments allow to extend the 
dynamic range of CMB observations to larger multipoles 
with the respect of WMAP7, thus measuring the damp- 
ing tail of the CMB angular power spectrum. While 
SPT probes the small scales in the range of multipoles 
650 < £ < 3000, the ACT telescope spans a range of mul- 
tipoles that goes up to £ = 10000, although the signal at 
i ^ 3000 is dominated by the power coming from extra- 
galactic point sources. For this reason, for ACT we only 
consider the less contaminated 148 GHz spectrum up to 
^max = 3300 to perform cosmological parameters extrac- 
tion. In order to account for the foreground contributions 
at £ < 3000, we add three extra amplitude parameters: 
the Sunyaev-Zel'dovich (SZ) amplitude Agz, the clus- 
tered point sources amplitude Ac and the amplitude of 
Poisson distributed point sources Ap. We consider for 
both ACT and SPT experiments a joint amplitude pa- 
rameter for each component and the templates provided 
by |6]. No SZ contribution is considered for WMAP7, 
as explained in the analysis performed by [H]. We have 
however verified that different choices for the foreground 
templates has negligible effect on the constraints of cos- 
mological parameters and produces minimal effects on 
foreground parameters. 

For what concerns the SDSS LRG7 catalog, we chose 
to consider data only in the linear scales regime, i.e.. 
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up to = 0.1 h Mpc-i. Indeed, HALOFIT, the CAMB 
routine that should correct for non-linearity effects at 
the smallest scales, is tested only for a smooth primor- 
dial power spectrum and therefore is not appropriate for 
dealing with power spectra with features, like those con- 
sidered in the present analysis. 

Regarding the prior on the model parameters, we im- 
pose flat priors on Ub, ujc, 0, r and Ug and a logarithmic 
prior on As- We check a posteriori that these priors result 
to be much wider than the corresponding posteriors and 
thus their upper and lower limits do not affect our final 
results. The priors on the step parameters need however 
to be discussed in more detail. The parameter b controls 
the position of the oscillations in fc— space. Larger val- 
ues of h correspond to "later" phase transitions and thus 
move the oscillations towards larger scales (smaller val- 
ues of k and t). Viceversa, smaller values of b shift the 
oscillations in the direction of large wave numbers. As a 
rule of thumb, we note that the peak in the oscillations 
is located at fc ~ 0.015 Mpc^^ {£ ~ 200) for b = 14.5, 
and that it is shifted down (up) by roughly a factor 2 
in k for each 0.1 increment (decrement) in b. Thus, out- 
side of a given range in 6, oscillations are moved to wave 
numbers that are not probed by observations of the CMB 
nor of large scale structures. Based on the considerations 
above, we initially choose a flat prior for b in the range 
13.5 < 6 < 15.5, that conservatively encompasses the 
whole range probed by the WMAP, ACT, SPT and LRG 
datasets. We use this prior for the CMB only dataset. 
Then, in view of the results of the first Monte Carlo run, 
we also consider a restricted prior 14.2 < & < 15, that we 
use for the analysis of the enlarged dataset. We have also 
explicitly checked that adding, in 6 = 13.5 or & = 15.5, 
a step-like feature with c — 1Q~^ (a value already large 
enough to produce, on average, oscillations that are at 
variance with observations [22]) and d = 3 x 10^^ (the 
median point of our prior) to the WMAP7 best-fit model 
produces no appreciable effect (at least within CAMB's nu- 
merical precision) in the CMB spectrum up to £ = 3000 
nor in the matter power spectrum between k = 0.02 and 
fc = 0.1 Mpc^^. For what concerns c and d, parameter- 
izing the height and width of the step respectively, we 
choose a logarithmic prior for both of them, i.e., a uni- 
form prior on log c and log d. The reason for this choice is 
that we want for these parameters to, potentially, assume 
values spanning several orders of magnitude with equal 
a priori probability. Indeed this is accomplished using a 
logarithmic prior that naturally assigns equal probability 
to each decade. In particular, we take —4 < logc < — 1 
and -2.5 < \ogd < -0.5. 

We derive our constraints from parallel chains gener- 
ated using the Metropohs-Hastings algorithm. We use 
the Gelman and Rubin R parameter to evaluate the con- 
vergence of the chains, demanding that R — 1 < 0.04. 
We note that there are some issues related to the fact 
that some of the posteriors do not vanish at infinity; we 
address them in section IIVI . 



IV. RESULTS AND DISCUSSION 

First of all, we check how the constraints on the six 
"vanilla" parameters are changed when the possibility of 
having features in the primordial power spectrum is con- 
sidered. To this purpose, we show in Tab. |l] the mean 
of the posterior distribution, and the corresponding 68% 
credible intervals, for each vanilla parameter, as well as 
for some derived parameters (most notably the reion- 
izaion redshift Zrc, the age of the Universe, the Hubble 
constant Hq), and compare them with the corresponding 
values found by the WMAP7 team in their analysis of 
the ACDM model; however, in order to have the max- 
imum homogeneity between datasets, we consider the 
dataset dubbed "WMAP7-hCMB" in the LAMBDA web- 
site that includes, in addition to WMAP, also the data 
from small scale CMB experiments. We note that the 
uncertainty on the determination of the vanilla parame- 
ters is not degraded when features are included (and it is 
actually better for the CMB+LRG dataset, although this 
should probably be ascribed to the inclusion of additional 
data). The mean values found for the features model are 
all within one sigma of the corresponding ACDM values, 
with the partial exception of the baryon density ajf,/i^, 
whose mean lies at the edge of the WMAP 68% credible 
interval. We argue that this is due at least in part to the 
fact that some of the primordial spectra considered here 
mimic the presence of a negative running in the spectral 
index (see below); the lower value of ujbh'^ is thus due to 
the correlation with the running. 

Now let us analyze the results on the primordial spec- 
trum parameters from one dataset at a time, starting 
from the CMB-only dataset. We recall that this analysis 
assumes the prior 13.5 < 6 < 15.5. We first focus the 
best-fit parameters, i.e., the parameter values that max- 
imize the likelihood. In the second column of Tab. Ullwe 
show the best-fit values, for this dataset, of the primor- 
dial spectrum parameters. In the best-fit model, the step 
in the primordial spectrum is located in 6 = 14.66, consis- 
tently with previous studies [20l [22] . We explicitly show 
the projected hkelihood, as well as Ax^ = X^ — Xmin' ^ 
function of b in the two panels of Fig. [T] It is interesting 
however that a distinct although low peak is present in 
the likelihood in & ~ 14. We also note that the does 
not become arbitrarily large (i.e., the likelihood does not 
asymptotically vanish) far from the minimum, but in- 
stead tends to a constant value. This is related to the fact 
that, as explained to the previous section, towards the ex- 
tremes of the b range the oscillations are moved outside 
the range of scales probed by the dataset, and thus the 
model becomes completely equivalent, as long as data 
fitting is concerned, to ACDM. Thus from the plot we 
can roughly estimate that the best-fit at 6 = 14.66 repre- 
sents a Ax^ ~ 7 improvement with respect to the ACDM 
best-fit, while for the model with b = 14, Ax^ ^ 3.5. The 
primordial power spectra correspoding to the two min- 
ima in the x^ are shown in Fig. [2j On the other hand, 
models with 14.1 ^ 6 ^ 14.5 perform worse with respect 
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TABLE I. Posterior mean for the vanilla cosmological parameters. The errors refer to 68% credible intervals. 



Parameter 
lOOfib/i^ 

100 6* 



Age [Gyr] 
Ho [km s~i Mpc-i] 



ACDM (WMAP7+CMB 



Features (CMB 



Features (CMB+LRG 



2.253 ± 0.054 
0.1103 ±0.0052 
1.0396 ±0.0025 
0.088lS°}| 
0.962lg;gl^ 
2.15 ±0.11 
13.72 ±0.12 
10.5 ± 1.2 
71.4 ± 2.4 



2.204 ±0.044 
0.1125 ±0.0050 
1.0409 ±0.0016 
0.086 ±0.014 
0.959 ±0.014 
2.18 ±0.08 
13.81 ±0.09 
10.5 ± 1.2 
69.9 ±2.3 



2.215 ±0.037 
0.1122 ±0.0029 
1.0414 ±0.0015 
0.087 ±0.015 
0.959 ±0.011 
2.19 ±0.07 
13.78 ±0.07 
10.5 ± 1.2 
70.3 ± 1.3 



^ Posterior mean for the parameters of the ACDM model, using the WMAP7+CMB dataset from |http:/ /lambda. gsfc.nasa.gov[ 
^ Posterior mean for the parameters of the features model, using the CMB dataset. 

Posterior mean for the parameters of the features model, using the CMB+LRG dataset 

ko = 0.05 Mpc~^. The ACDM value reported in the Lambda website is for kg = 0.002 Mpc~^ but has been rescaled to ko = 0.05 to 
allow for comparison. 




FIG. 5. One-dimensional posterior probability density for the step parameters from the CMB dataset. 




TABLE II. Bestfit values. 



Parameter 


CMB 


CMB+LRG 


b 


14.66 


14.66 


logc 


-2.65 


-2.85 


logd 


-1.42 


-1.57 


Us 


0.946 


0.958 


ln[10^°A.] 


3.08 


3.06 


-log(£) 


3765.8 


4043.7 



to ACDM. 

For a better understanding, we also show in Figs. [3] 
and |4] the comparison between the WMAP7 best-fit, the 
two models with features corresponding to the two peaks 
in the likelihood seen in Fig. [l} and the data present in 
the CMB dataset. It is clear from these plots (especially 
from the plot of residuals shown in the lower panel) that 
the model with b — 14.66 improves over ACDM by be- 
ing able to fit the two outliers in £ = 22 and £ = 40, 
thus confirming our previous findings |22| . The inter- 
pretation of the peak in 6 = 14 is more puzzling; by 
looking at the lower panel of Fig. [Sj however, it can 
be seen that the CMB spectrum for this model resem- 
bles what it would be obtained by adding a negative 
running dng/dlnk to the scalar spectral index. Thus 




FIG. 6. Joint two-dimensional posterior for (6, logc) using 
the CMB dataset. 



this result could be reminiscent of the WMAP7 prefer- 
ence for a negative running, that is indeed even more 
pronounced when high-^ data are added to the analy- 
sis [T]: dris/dlnk = -0.034 ± 0.026 (WMAP7 only) and 
dnjdlnk = -OMltoml (WMAP7+ACBAR+QUaD). 

We now turn to the posterior distributions. In Fig. 
[5] we show the one-dimensional posteriors for 6, log c 
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and logd. The posterior for b still shows the two peaks 
in & ~ 14.7 and 6 ~ 14 that were present in the likeli- 
hood. The largest value at the edges of the prior range 
is due to a volume effect, since the one-dimensional pos- 
terior is obtained by marginalization (as opposed to the 
one-dimensional likelihood that was obtained by maxi- 
mization). On the other hand, the probability density 
for 14.1 < 6 < 14.4 is practically equal to zero. For what 
concerns logc, as it could be expected, "large" values 
are disfavored by the data (as they produce large - in 
amplitude - oscillations that cannot, on average, be rec- 
onciled with observations) while for smaller values the 
posterior tends to a constant value as the oscillations 
become so weak as to be practically undetectable for 
the current experimental precision and thus the value 
of c becomes unimportant. As already noted, a posterior 
with this characteristic that extends, in principle, down 
to logc = — oo, cannot be properly normalized (since the 
corresponding probability mass is infinite) and, as a con- 
sequence, credible intervals are ill-defined. One could be 
tempted to impose a lower cut-off but then the credible 
intervals will end up depending on the choice of the cut- 
off itself, so this should be avoided, at least in the absence 
of a clear physical reason for doing so. 

We can still, however, compare probability densities, as 
well as probabilities integrated over finite intervals, since 
probability ratios do not depend on the overall normal- 
ization. We can use, as a benchmark value to compare 
the constraining power of different datasets, for example, 
the value of log c where P(log c) is half of its asymptotic 
value for log c — > oo. This should not be taken as an "up- 
per limit" in the common sense of the word, but as said 
is a useful tool for comparison. In the case under consid- 
eration, we estimate that this happens for logc — —2.32, 
or c = 4.8 X 10"3. 

For comparison, the corresponding value that we had 
previously found using WMAP7 and ACT data only was 
logc = —2 (33. We also show, in Fig. [6] the two- 
dimensional posterior P{b, logc) where it is clear that 
probability is concentrated in two distinct, disconnected 
regions. One corresponds to models with b ~ 14.7 and 
logc ~ —3, while the other to models with 6 ~ 14 and 
logc located more towards the edge of the prior range, 
logc ^ —3.5. Finally, we examine the posterior for logd. 
This is in part similar to the posterior for log c, once one 
recalls that small values of log d produce a steep step in 
the potential and consequently large oscillations, so one 
should expect the probability to go to zero for small val- 
ues of log d, as it is. However, in this case, the posterior 
range is not wide enough to see the asymptotic part, for 
logc? — > 00 (where ACDM should be recovered), of the 
distribution. 

The fact that the posterior is bimodal in b creates some 
difficulty for the Monte Carlo, as the chains cannot eas- 
ily jump from one peak to the other, and thus take a 
longer time to sample satisfactorily the actual distribu- 
tion. For this reason, in our second Monte Carlo run, 
using the CMB-I-LRG dataset, we have decided to con- 



centrate on the region of the peak at 6 = 14.66 and im- 
pose the prior 14.2 < 6 < 15. We find that the best-fit 
for this dataset, shown in the third column of Tab. |lT] 
has still b = 14.66. In Figs. [7] and [8] we compare the 
one-dimensional likelihoods and posteriors, respectively, 
for the step parameters in the CMB-I-LRG dataset with 
those obtained previously with the CMB dataset. In or- 
der to allow for comparison, the distributions for the 
latter have been obtained by imposing a posteriori the 
condition b > 14.2 (which, in practical terms, that we 
have discarded all samples with b < 14.2, and reanalyzed 
these new chains from scratch). We find that there is 
practically no difference with respect to the position of 
the oscillations (which makes sense, since this is driven 
by the requirement of fitting the outliers in the WMAP7 
data at relatively low ts). The amplitude of the oscil- 
lations is slightly more constrained, with the posterior 
going down at half of its plateau value at logc = —2.48 
(c ~ 3 X 10^'^). The posterior for logd is also slightly dif- 
ferent, as it shows a more distinct peak in correspondence 
of the best-fit value logd ~ —1.5. 



V. CONCLUDING REMARKS 

We have studied cosmological models with a step-like 
feature in the inflationary potential. Such a feature 
would produce oscillations in the primordial spectrum 
of scalar perturbations, whose presence can be tested 
through the analysis of CMB and large-scale structures 
data. We have found, consistently with previous studies, 
that in these models the agreement with the CMB data 
is improved, with respect to the ACDM model, when the 
oscillations are placed in such a way as to match the two 
outliers in the WMAP7 spectrum at ^ = 22 and ^ = 40 
(in particular, the changes by Ax^ — 7). The poste- 
rior probability also has a maximum close to this point, 
corresponding to 6 = 14.66, while it clearly shows that 
oscillations in the range 14.1 < & < 14.5 are currently 
forbidden by the data. The possibility of no oscillation 
at all is still, however, perfectly consistent with the data. 
In conclusion, although multifield inflationary models can 
definitely reproduce the two glitches in the WMAP7 tem- 
perature spectrum, current data are not yet constraining 
enough to allow to discriminate between these models 
and the standard inflationary scenario. 
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FIG. 7. Comparison of the model likelihood as a function of the step parameters, obtained from the the CMB (red solid line) 
and CMB+LRG (blue dashed) datasets. 




FIG. 8. Comparison of the one-dimensional posterior probability density for the step parameters from the CMB (red solid line) 
and CMB+LRG (blue dashed line) datasets. 



